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ABSTRACT 

Time evolution of the decay process of unstable particles is investigated in field 
theory models. We first formulate how to renormalize the non-decay amplitude 

P-i. 

D , beyond perturbation theory and then discuss short-time behavior of very long-lived 



particles. Two different formalisms, one that does and one that does not, assume 
existence of the asymptotic field of unstable particles are considered. The non-decay 
amplitude is then calculated by introducing a finite time resolution of measurement, 
which makes it possible to discuss both renormalizable and non-renormalizable decay 
interaction including the nucleon decay. In ordinary circumstances the onset of the 
exponential decay law starts at times as early as at roughly the resolution time, but 
with an enhanced amplitude which may be measurable. It is confirmed that the 
short-time formula 1 — Ft of the exponential decay law may be used to set limits 
on the nucleon decay rate in underground experiments. On the other hand, an 
exceptional example of S-wave decay of very small Q-value is found, which does not 
have the exponential period at all. 
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1. Introduction 

It is an old problem since the classic works of Dirac, Weisskopf and Wigner how 
to describe the decay process of unstable state, following the principles of quantum 
mechanics. It has been repeatedly questioned how the exponential decay law is mod- 
ified both at late and early times. New motivation for investigating this age-old 
problem has been added when relevance of the non-observation of the nucleon decay 
is addressed 0, ||^ and whenever high resolution experiments with new technol- 
ogy are performed to search for deviation from the exponential law]^], [0. A 
conceptual problem, called the quantum Zeno paradox]^, also renewed interest in 
this problem, the central issue here being whether the initial decrease of non-decay 
probability proceeds via the square of time, t^, which appears inevitable in quantum 
mechanics, but not in accord with the exponential law. 

A field theoretical description of unstable particle decay has complication due to 
renormalization. It was pointed out in ref that the usual mass, coupling and wave 
function renormalization is not sufficient to remove all infinities that appear in the 
non-decay amplitude at finite times, if one works in perturbation theory. We shall 
discuss below how this problem may be resolved. 

In the present work we first discuss an exact formulation of time evolution of 
unstable particle decay without assuming existence of the asymptotic field of unsta- 
ble particle. We then show how renormalization is performed beyond perturbation 
theory, both in this formalism and in a more conventional approach that assumes 
existence of the asymptotic field. 

We next investigate how time evolution of the decay proceeds when viewed with 
finite time resolution. Finite time resolution can be regarded as an effective means of 
introducing a resolution in actual observation. In all applications of practical interest 
the resolution time At is limited from both below and above, ^ > At > where 
M is a physical energy scale of the process such as the mass of the decaying particle 
and r is the decay rate. If the time scale one wants to probe into is larger than 
Y, description of the decay as a function of time effectively loses its meaning, and 
one looks into some other measurable quantity such as the energy peak of produced 
decay products directly. If the time scale is smaller than -p, the concept of particle 
picture is doubtful. It is thus sufficient to be able to compute the time evolution 
of the non-decay amplitude only for a coarse grained resolution in the time range 
above. 
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The decay law then differs in fine detail, depending on the amount of resolution. 
In ordinary circumstances the exponential law becomes excellent at times as early as 
of order At. The absolute value of the non-decay amplitude is however enhanced by 
a resolution dependent amount, ~ 1 + | TAt for < At ^ ^ , and 1 + for 
At ^ 2^, where Q is the Q value of the threshold. The short-time behavior of the 
probability at t < At shows a typical time dependence of 1 — t^ , characteristic 
of the quantum mechanical law. When At becomes large and approaches the lifetime 
Y, time evolution shows a large deviation from the exponential law. 

An interesting exception on observability of deviation from the exponential law 
occurs when Q <r, which becomes possible in the Q ^ limit of the S-wave decay; 
r ^ (coupling) X \JMQ . The exponential period may not exist at all in this case. 
Whether this case is realized in nature is yet to be found. 

Another common view of the unstable particle decay is to consider as inseparable 
the formation process of unstable state created from stable many particle states. [§] 
This view introduces a new time scale of formation time interval. But if the time 
resolution At we consider here is much larger than this formation time interval, then 
effects of finite time resolution dominate over the formation effect. This inequality 
of two time scales usually holds in most practical cases. This is fortunate because it 
is usually difficult to know the formation time interval. 

A field theoretical approach discussed in detail in Appendix makes it possible to 
treat the region of theoretical interest, t -C 1/M. The infinitely good time resolution 
is realized in this formalism; At = formally. It is found that after renormalization 
of the proper self-energy the non-decay probability has some peculiar behavior in 
the t — >■ limit, at least in perturbation theory. Moreover, the limiting behavior is 
not universal. In a renormalizable case examined there is a linear rise oc 1/t in the 
time range, ^^^^^ > t ^ ^~ '^'^^ ^ 17i ■ ^^^^ other hand, one has the linear decrease 
of 1 — 0[t] in a super-renormalizable model. These results cannot be extrapolated 
to the true t ^ limit, because one must go into a non-perturbative region. 

It thus seems that one cannot answer the quantum Zeno paradox of unstable 
particle decay in field theory that assumes the asymptotic field. The true t ^ 
limit, or At ^ limit, is difficult to deal with in perturbation theory, the only 
scheme that one can compute to a high precision, because in approaching t — » one 
has to know the behavior of Green's function in ever increasing high energy. 
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2. Projector formalism and decay law in quantum mechanics 

We shall use an exact integro-differential equation for the non-decay amplitude, 
(l|e~*''^*|l) , with |1) the inital state of unstable particle. Following Peres0, we 
introduce the projector onto the initial unstable state, P = |1)(1| , and decompose 
the total Hamiltonian H as 

Ho = PHP + (1 - P)H{1 -P), V = H-Ho. (1) 

Define energy eigenstates \a) of Hq; Ho\a) = Ea\a) . It is then easy to confirm that 
the decay interaction V operates only between the initial state |1) and its orthogonal 
compliment \n); Vi„ = V*-^ ^ , Vmn = Vn = . We use the intermediate Roman 
letters such as m , n to denote eigenstates projected by 1 — P. The interaction V is 
defined here in a novel way; it depends on the prepared initial state. Although this 
looks odd, it makes possible some exact analysis, as will be made clear shortly. 

We work in the interaction picture and expand the state at a finite time t, using 
the basis of the eigenstate of Hq ; \ip)j = J2a (^a(t) \a) . We thus write the time 
evolution equation for the coefficient Ca(t); 

tCm = Knl e^(^™-^l)* Ci , 2Ci = ^ Vi„ e-^^^""^^)* Cn . (2) 

n 

The non-decay amplitude is related to this coefficient by 

(l|e~*^*|l) = e-'^^' (l|e*^°V*^*|l) = e~'^'' = e-'^''ci{t) . (3) 

A closed form of equation for the non-decay amplitude Ci{t) = c{t) then follows; 

c{t) = - r dt' (3{t - t') c{t') , (4) 
Jo 

roc 

(3{t-t') = {l\Vi{t)Vj{t')\l) = duja{u)e~'^^-^^^^'-'\ (5) 

J Uc 

(yi^^) = ^ ^i^^ - Era)\Vimf' . (6) 
m 

Here Vi{t) = e^^°^ y ^-^Hot jg ^j^g decay interaction written in the interaction picture 
and the spectral function cr{uj) characterizes interaction between the unstable particle 
and decay product particles. The initial condition Cm(0) = was used to derive the 
equation for c(t), and uJc is the threshold for the state \n). 

The familiar golden rule of perturbation theory immediately follows if one ap- 
proximates c{t') in the right hand side of eq.(^ by unity and take the infinite time 
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limit; |c(t)p ^ 1 — 27Ta{Ei)t . Thus the perturbative decay rate F is 2TTa{Ei). We 
would however like to elucidate the time evolution in finer detail, not relying on 
perturbation theory in this manner. 

The standard technique to solve this type of integro-differential equation is to 
use the Laplace transform, and one finally obtains the non-decay amplitude in the 
form, 

l roo ^-i{uj-Ei)t 

'^^^ = 2m J-oo'^"' -iu + E,-Giuj + tO+) 

1 POO 

= / dcjF(cj + zO+)e-^("-^i)*, (8) 

2iri J-oo 

G{uj + tO+) = r dJ , ""^""'^^^ = I r rfr/?(r)e^^('^-^^)-°''^ . (9) 

Ju)c. UJ' — UJ — 20+ JO 

The initial condition c(0) = 1 was imposed in this derivation. 
The analytic property of the function, 

^''-'" -.+E!-G(.r 

is evident; F{z) is analytic except on the real axis with the branch cut starting from 
a threshold, z > uoc- This makes it possible to write the amplitude in a convenient 
form, 

/■oo 

c(t) = / dcjiJHe-*("-^^)*, (11) 



^ ' 27rz (cu - El + n(cu))2 + (7r(TH)2 ' ^ ' 

where n(ti;) was defined by the boundary value of G{z)^ when z approaches the real 
axis, G{uj + zO"^) = n(ci;) + ma{ijj) . 

It is well known that in the cut uo plane there is a pole in the second Riemann 
sheet, near (and below) the real axis if the decay interaction is weak enough. The 
pole location z is determined by 

z-E^ + Gn{z)=Q, (13) 

where the analytic function G{z), hence F{z), is extended into the second sheet by 
Gii{uj — zO^) = Gi{uj + i^'^) through the branch cut at > uj,.. Dominance of 
this pole term is equivalent to taking a Breit-Wigner form for H{uj) in eq.(0) and 
neglecting the threshold value uo^ by the replacement, ujc ^ — oo . The equation 



5 



([T3|) gives perturbatively an imaginary part of the pole, —iT/2 = — Z7ra(£'i) . The 
ideal exponential decay law |c(t)p oc e~'"* then follows, but this approximation fails 
both at early and late times. 

At late times the threshold behavior of (t(c<j) taken to be cr{uj) = c{uj — uJc)°' 
dictates the power law |c(i(:)p oc t"^*^""*"^-' [jlO|- The transition time tp from the 
exponential to the power decay may be estimated by equating the amplitude in the 
power period, 

c(t)power - - e-(-^^ , (14) 



to the exponential e '"'/^ to give |TT]| 

*^ " f ( cr(afi)r^.j • <^^) 

Here T{z) is the Euler's gamma function, and Q = Ei — uJc is the Q-value for the 
decay. 

By using the "elastic" unitarity relation, 

F{uj + iO+) - F{tu - iO+) = 2wH \F{uj + zO+)p , (16) 

equivalent to eq.(|l2D, one may deform the contour of integration on the real axis 
into a sum of two contours, the one around the pole as shown in Co in Fig.l and the 
other along Ci; 

c(t) = —([ +f ) dz F(z)e-'^'-^'^K (17) 
2m \Jco JCi J 

The limiting power law formula, eq.(|n|), is a result of an approximation to the 



contour integral along Ci, which can be used in the large time limit. In Fig. 2 we 
compared eq.(p!4D and an exact numerical computation of the contour integral along 
Ci, which shows that the formula ([T^) is an excellent approximation. 

As an application we used this formula to estimate the transition time to the 
power law decay for the tt /xz/^ decay. We use 



The result is 



a(u;) = ^{—f ^^ "^V , (18) 
^ ^ 27r^mJ (1-4)2' ^ ' 

^^ = lGlf^m.mlil-^r. (19) 



2 / TT , , 4 , "^u ^ 9 / X 9 \ 280 , , 
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Since e~'"'"*p is too small, it is difficult to observe the power law decay in this example. 

On the other hand, at early times the high frequency behavior of the spectral 
function cr{uj) becomes important, on which we shall mainly focuss. One naively 
expects that the short-time behavior exhibits a deviation from the exponential law 
in the form of 

|(l|e-'^*|l)P ^ 1 -t' [{1\H^\1) - {l\H\lf) . (21) 

Thus, quantum mechanics appears to predict the quadratic form of deviation in the 
t —>■ limit. More precisely, our solution (|TTp gives 

|c(t)|2 ^ l-t2A(cj-^i)2, AA2= / dujE{uj)A^-[ \ dujH{uj)A] .(22) 

As will be made more explicit shortly, field theory models, worked in perturbation, 
give at best (in the super-renormalizable case) the high energy behavior of 
const /w^, as u; oo. The coefficient of order term is then infinite in field 
theory. It is however not clear whether the power series expansion inside the lo 
integral for estimating t — > limit is mathematically permissible or not. 

3. Field theory models and renormalization 

The projector formalism thus presented recasts the non-decay amplitude into the 
form known in an exactly solvable field theory model, the Lee model. We briefly 
summarize this connection below. For general renormalizable field theory models 
there is however a great difference from the Lee model, in that the spectral weight 
oiiS) in the projector formalism cannot be arbitrarily given and is fixed by the theory, 
which however is difficult to extract beyond perturbation. We shall nevertheless use 
the form of renormalization suggested from the study of the Lee model for unstable 
particle decay. [ 13 1 Renormalization in a more conventional approach that assumes 
existence of the asymptotic field of unstable particle is dealt with in Appendix. 

Let us ffist note that the projector formalism may be recast into an equivalent 
quantum mechanical model. We imagine that the decay process is a transition 
between the states, |1) . All other states that decouple from the decay 

process may be ignored, which should not pose any problem for our purpose. With 
Ex > uJci the Hamiltonian h of this model is of the form, 

/•oo /"OO I 

h = Ei\l){l\+ duuj\uj){uj\+ duj^a{uj){\uj){l\ + \l){uj\) . (23) 

JUJr JUJr 
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We shall show the equivalence after we introduce a field theory model. 

We assume that the states relevant to the decay process are given by a Fock space 

of 

\l)=a^O), \uj) = b^{uj)\0) . (24) 

We may then extend the above model to an equivalent quantum field theory model 
(EQFM). The model is defined by 

H = {Ei + 5E)a)a+ duo uoh\uo)h{uj) + —= duo ^ a{uo) [ah\uo) + a)h{uo)) . 

(25) 

Here 5E and Z are renormalization constants needed to compensate against a bad 
high frequency behavior of (j{uo). 

This model has an exact conservation law; the sum of the total a— and 6— particle 
numbers is conserved. This enables one to solve this model in each sector of a given 
conserved particle number, separately. The consequence of the projector formalism 
follows by identifying the projector P = a"''|0)(0|a in this field theory model. 

The freedom of choosing renormalization constants, Z , 6E, are fixed later by a 
renormalization prescription. Allowing for their existence is crucial to renormaliz- 
ability. If Z is complex, one may redefine the phase associated with the state |1) 
such that Z may be taken to be real and positive. Although written in a different 
form, it should be evident that this model is in essence equivalent to the Lee model 
for the case of the vanishing particle mass, rriN = 0. In the Lee model for unstable 
particle decay this factor Z is infinite in the limit of local field theory (when the so 
called cutoff function is unity) . |T3| 

We assume the canonical commutation relation; 

[b{uo),b\uo')] = 5{uo~uo'), [a,a^] = l. (26) 

Since the continuous states \uo) are two-body states of the decay product, assumption 
of the canonical commutator for b{uo) is highly non-trivial. It is conceivable that 
there is some correction to this relation. Neglect of such possible correction should 
be regarded as an essential part of the approximation to the full theory. Furthermore, 
there are parts of the full-fledged Hamiltonian such as a b^{uo) that are not included 
here, because they are not relevant to one-particle decay process. When one discusses 
many-body problem of unstable particles, they may become relevant. 
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We also note that in our EQFM approach we do not assume existence of asymp- 
totic fields. This is welcome because the asymptotic field corresponding to unstable 
particle is not well defined in a strict sense. 

To show that time evolution of the state |1) in EQFM is equivalent to the one in 

the projector formahsm, we first seek energy eigenstate Ici;)^ with {H — ui) \ijj)s = , 
in the form of scattering state; we find that 

\uj)s^B\ujm, (27) 

/ I o) 1 \/ C!'(uj)cr(u') ^ \ 

BtM = 6tM + F(. + «-) + ,.'1J_LL1,V)) . (28) 

F{z)-^ ^-z + Ei + SE-^ r cLu . (29) 

There exists an important relation of Hamiltonian equivalence, 

r duuB\u)B{u)). (30) 

Furthermore, two descriptions, either in terms of (a , 6(a;)) or (5 (a;)), are related by 
a canonical transformation, satisfying 

{B{uj),B\J)\ = b{uj-J). (31) 

We further note a remarkable operator inversion; 

a^^-^ / (iu;'\/(7(a;)F*(a; + iO+)S'f(u;) , (32) 

1 1-^ . \AM^(^') i^*(c^' + iO+) , , 

The proof of all these is based on the analytic property of F[z) having the cut 
on the real axis from z > uJc, with the discontinuity formula of 

2-ni 

F(u; + «0+) -F(u;-^0+) = — cT(a;) |F(a; + iO+)|^ for a; > a;c . (34) 

Zj 

Its asymptotic behavior, F{z) — \ i is also important. This asymptotic behavior 
holds only when a finiteness condition for 

E,^8E-\ Hdu;^, (35) 

cu , , cr(uj') , , 

/ <i-^ ,/ , , 36 

Z Lc -UJ- iO+) ' ^ ' 
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is fulfilled. 

Time evolution of the state |1) is derived, by noting the completeness ||I4 



du\uj)ss{uj\ = \l){l\+ du\uj){uj\. (37) 



We find from 



1)= / diUs{co\l)\uj)s, (3^ 



that 



1 



e-*^*|l) = -— / dujJa{uj)F*{uj + iO+)e-'^'\uj)s, (39) 



(l|e~'^*|l) = (0|ae-'^*at|0) = ^ I" doo a{uj) \F{uj + iO+)\^ e'''^' . (40) 



Z 

Due to the cut structure of F{z), one can write this as 

1 

(l|e-'^*|l) = — / rfcue-'"*F(cu + zO+) (41) 
zm j~oQ 



\D{uj)\^ = (uj-E,-6E + ^V r dcu'^^] 

Y Z JuJc ijj' — ijj j Z^ 



\2 



.(43) 



This result coincides with the result of the projector formalism (|TTp except for the 
Z factor and 6E. 

There are various ways of how to specify the renormalization condition. We adopt 
an interesting proposal using the residue of the physical pole of unstable particle in 
ref |T^]; |-F(2;)| ~ ■> where E^, = E,. — iEi is the pole location in the second 

sheet, and is some finite constant. The renormalization condition then reads as 

Z Ju, ^ ' \uj-E,\^ I Z^ \Ju., ^ ' \uj-E,\^ I ^ ^ 

For a given high frequency behavior of <y{oj) renormalizability demands that all 
infinities, if there is any, must be eliminated by two infinite constants, 5E and 
Z, considering that the renormalized and observable quantity E^, is finite. The 
quantity Z \D{ijj)\'^ appearing in the non-decay amplitude (^2]) should then be finite 
in renormalizable theories, possibly with an exception of an infinite phase factor that 
does not contribute to the non-decay probability. The renormalization condition 
may be called the on-shell renormalization condition. 
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When the integral 

A,^r,^^ (45) 

J UJc ^ 

is finite, the infinity appears only in the mass shift, 5E. On the other hand, when 
this integral is infinite, one finds that 

N (uj - E - E"^ 

dua{u) ^^ ^V.M • (46) 



1 — N Jujc \lj — E^^ 

The criterion for existence of the remaining infinity Z is that the integral Ai is 
infinite. We note that the on-shell renormalization condition (HI) is necessarily non- 



perturbative, in the sense that the solution Z does not behave like Z = 1 + 0[a]x 
an infinity, as might have been expected. If N ^ 1 + 0[a] , it behaves rather like 
Z = O [cr°] X an infinity. 

We now show some examples of the spectral function (j{uj) in perturbation. Three 
models are taken; a superrenormalizable model of boson decay defined by an interac- 
tion Lagrangian density, Cb = ^ cpx^ , a renormalizable model of fermion-pair decay, 
Cf = g ipipip , and an effective model of proton decay, Cp = e(a + 675)^7? + (h.c.) , 

2 2 

where a , b are constants of order, ^^^p . Here mx is the mass (~ 10^^ GeV) of 
heavy X boson mediating baryon non-conservation. Let us first note that in the 
projector formalism, or more simply in EQFM, there is no Z-diagram as shown in 
Fig. 3, contributing, and one has to consider the direct diagram alone. Perturbative 
computation in the couphng, ju and g, then gives the spectral function a^uj) to lowest 
order; 



Here M and m denote the parent mass and the daughter mass, respectively. (In the 
case of the proton decay we can take the electron mass vanishing such that m is the 
pion mass.) The lowest threshold uJc = 2m (or = m for the proton decay). 

We point out that in the case of fermion-pair decay there is a quadratic divergence 
for the self-energy diagram of — > tptp —>■ if , hence the spectral function increases 
as (Tfiuj) oc cu^ for u ^ 00 . This bad high energy behavior in the renormalizable 
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model introduces infinities in the integral, du^j^ for G{z), proportional to 
z° , z , z"^ , and the term independent of z is eliminated by the mass renormalization. 
The remaining integral contains an infinite Ai term, and the renormalizability of this 
particular field theory model is determined whether 

Ipfrf,'-;^ (50) 

is finite or not. The fermion-pair decay model here indeed obeys this finiteness 
condition. 

In Appendix we show that if one takes the initial one-particle state at rest given 

by 

\1) = -i J d X fo{x,ti) do ^{x,ti)\0) , foix,ti) = -j= , (51) 

where if{x) is the Heisenberg operator of unstable particle (whose existence is as- 
sumed) and |0) is the true vacuum state, then the non-decay amplitude is related to 
the ordinary Green's function; 

^^|g-ii/(t-tO|i) = (!^L±ill! f d?x [ d^Xi{0\T{ip{xMx,))\0) . (52) 



2M 

Assuming a convergence of integral, one can write this as 

r<ic.^^ A(.V,0-). (53) 

where A is the usual Green's function expressed in terms of the Lehmann spectral 
function pL^cr^); 

A{u') = -z / ci^xe^-*(0|T(^(a;)^(0))|0) = T da' ^'^f},^^ - (54) 

It should be evident that the non-decay amplitude is finite after renormalization if 
the Green's function is renormalizable. This is in sharp contrast to the argument 
based on perturbative expansion 0. We shall discuss the problem of renormalization 
in the field theoretical framework further in Appendix. It is also true that the 
disconnected diagram is factorized in the non-decay amplitude, and one may ignore 
the disconnected diagram in subsequent computation. 

4. Field theory with finite time resolution 

To search for deviation from the exponenital decay at short times, one should 
keep in mind that following time evolution with infinitely good time resolution may 
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not be meaningful and certainly is not practical with a limited accuracy of time 
measurement. With this in mind we consider how the decay proceeds with a finite 
time resolution. As mentioned in Introduction, we take the resolution in the range, 
Y>At>j^. With this coarse graining one needs the spectral function a{uj) only 
in a relatively low frequency range, and the high frequency part is effectively cut off 
by a finite time resolution. This makes it possible to discuss the decay law even for 
non-renormalizable theories such as the four-Fermi type theory given by X boson 
mediated interaction. 

We time average the basic function — t') using a resolution function 5{t ; At), 

^(t; At) f« for |t| > At; ( 5(t ; At) ) =0, H dtd{t ; At) = 1 . 

\dt J J-oo 

(55) 

Examples of the resolution function are the step function and the Gaussian function. 

At _ Ul \ 1 ^2 

= S<^it;M)^ e-»^. (56) 



Another useful resolution function is the inverse Fourier transform of a Lorentzian 
form, 

giving in the real time 

^^(t;At) = i^(l + ^)e-W/^^ (58) 
With a finite resolution, the averaged function is 



/CO 
dr' (3(t')5(t' -r;At) 
-oo 



dwe-'(^-^^^^a{cu). (59) 
By defining the Fourier transform of the resolution function, 

/oo 
dte~'''H{t;At), (60) 
-oo 



one finds that 



Ty-^ l " ^ ' ( M M'^^l-jSi^At) 
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This formula has a simple interpretation; with a finite time resolution the original 
spectral function <j{uj) is modified to (t{uj) 6^-{uj — Ei ; At) . Thus, the direct diagram 
has a spectral function peaked at Ei of width 1/(2 At), which is cut off at the 
threshold ujc- 

Instead of j3{t — t') = {l\Vj(t)Vj{t')\l) , one may introduce the finite resolu- 
tion both at V/(t) and at V/(t') using the Fourier transformed resolution function 
Sf{u} ; At). In this case the Fourier transform of the resolution function for f3{t — t') 
is {Spioj ; At))^ ; the resolution function is modified, but it is still of the form of the 
resolution for f3{t — t'). Thus one may view our averaging as a coarse graining of 
local operators in quantum field theory. 

5. Application to unstable particle decay 

A universal effect of the resolution 6f is to cut off the high frequency contribution 
in the region of | a; — M | > . Furthermore, if At > ^ , the contribution near 
the threshold u ^ Q is also cut off. The effect of finite resolution then prolongs 
the epoch of the exponential decay, since the threshold behavior responsible at late 
times becomes less important. We shall estimate how deviation from the exponential 
law occurs at early times, assuming < At ^ ^ (decay hfetime). Obviously, for 
At > ^ the decay law is complicated and does not obey the simple exponential law; 
in this case one has to make indivisual analysis taking into account each specific 
situation. 

We first compute the amplitude of the exponential decay when this law holds. 
The magnitude of the exponential decay is determined by the residue of the pole in 
the second Riemann sheet. For this purpose we may use an approximate formula for 
the real part of the self-energy; 



where r is the decay rate in perturbation theory. We shall take a simplified form of 
resolution for analytical estimate. 




(62) 



M-;At) = .(^ 



x\). 



(63) 



For small g the pole term is then approximately given by 



F{z) 



(1 



1 - - r At ) 



1 



(64) 



z-M + 
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We have ignored a small imaginary part of the residue, which is of order V/M. The 
residue factor it! = ( 1 - | TAt is replaced by = ( 1 - for /^t < i^. 

The exponential period is approximately described by the probability function, 

\c{t)\'' ^ e-^\ = Max(At,-^ + =^). (65) 

The enhanced magnitude \ — becomes maximal around At = 2^ , with a value 
around Although time variation of the enhanced exponential decay is difficult 
to distinguish from the ordinary exponential law, it may be possible to observe the 
magnitude of the enhanced non-decay probability if one can measure the onset time 
of decay. 

Let us estimate the onset time of the exponential period, first by examining 
the non-decay probability at early times. Short time behavior of the non-decay 
probability is most sensitive to the high frequency behavior of (t{uj). By expanding 
g-i(a)-M)t^ one obtains for a small t 



poo / roo 

/ dcu Hr{u) {u-Mf -[ / duo Hr{u) {u - M) 

J2m \J2m 



(66) 



rr / (7(0;) 5f{u ■ At) 

where 11^ is determined by the dispersion relation from a{Lj)6F{uJ ; At). As a crude 
estimate one may take a truncated Breit-Wigner form. 



27r (u - M)2 + 



which can be used only when the Q value (— M — 2m) is not too small; 
2(M -2m)At>l. With FAt < 1, 

Thus, the coefficient of order t^ term formally diverges linearly as At — > . Since the 
power series expansion in time t is valid only for t < At, this does no harm practically. 
This simplified short-time formula is reasonably good, but the approxiamtion is not 
excellent, presumably reflecting dependence on a precise form of resolution function 
taken. 

More quantitative questions can be raised as to how the decay law is changed if 
one varies the resolution within physics mass scale. At > (M being the parent 
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mass). We have numerically computed the non-decay probability for the boson-pair 
and fermion-pair decay models, assuming various time resolution At. Examples of 
the result of this numerical computation for the case of the fermion-pair decay are 
presented in Fig. 4 and Fig.5. As qualitatively discussed above, deviation from the 
exponential law appears at short times of order At or less. At the same time the 
exponential term in intermediate times is enhanced by At dependent amount. In 
Fig. 4 we show the non-decay probability in the entire time range including three 
different time stages; early 1 — 0[t^] in the inlet of the figure, intermediate exponen- 
tial, and late power stages. In Fig.5 we show the short-time behavior in comparison 
to the limiting behavior; 1 — Olt"^] region using eq.(^), and the exponential region 
using a numerically computed pole location and its residue. In all these numerical 
computations we used the Lorentzian form of time resolution, eq.(|^. 

The transition time to the exponential period is crudely estimated by equating 
this short-time formula to the enhanced exponential formula; ^ — ^^ = -R^(l — Ft*) . 
It gives, with R'^ ^ 1 + ^ TAt , 

^ (tt - Vvr^ - 8) At ~ 1.8 At . (70) 

In our numerical computation this crossover of the two formulas of the regions does 
not actually occur. Nevertheless, the transition to the exponential period is given 
by a factor of a few times At, as seen from Fig.5. 

In estimating the short time limit of the non-decay amplitude, one needed the 
high energy behavior of the spectral function a{uj). Let us consider a cutoff theory 
of a finite time resolution At in which the Ai integral diverges in the At — > limit. 
If the time resolution At is small but is still much larger than the inverse of the 
cutoff scale 1/A of physical relevance, one may discuss the short-time behavior in 
the range t > ^. This consideration may be applied to non-renormalizable models of 
four-Fermi type, where A is the weak or the X boson mass. The use of perturbative 
formula for the spectral function (j{uj) limits the allowed range of the resolution At, 
since the integral Ai must remain of order unity. We shall use this constraint to 
examine the short-time limit of t ~ At, at which one may expect a deviation from 
the exponential law. 

For definiteness, we considered four examples; the fi decay, /i — ^ z/^ ez/g , the 
neutron decay, n — > peue , the vr decay, it Wi , and the proton decay, p e+7r° . 
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In all these cases we find, assuming /S.t > ^ and ^ 2^7' ^ universal result; 

may be approximated by „ ^A.^ in the integrand. Hence, the condition < 1 , 



smce 



gives 



Thus, the exponential law is good as early as the lifetime times of order, {■^Y , 
which is very small in all cases. 

Let us examine more closely the case of a small Q value; the neutron decay for 
At < 2^ and a hypothetical case of tt — > iDi decay, in which the mass difference is 
very small; m-,^ — nii ^ m^r . These are special examples of the small Q limit. First, 
let us discuss the tt decay. With Q = m^r — rrii and for the case of At < ^ , one can 
work out the constraint Ai < 0[1], to get 

pv...) ^ _ ^ (74) 

Because a typical dimensionless quantity in the left hand side is G\m\ fa 10~^^, 
this constraint is readily obeyed for At > l/rriT^. 

In the case of the neutron decay for a small Ai, one has, assuming At < 

Gl. i-mn-mp+l/2M , fm„+l/2At _ ^ _ jji\2 

^ / dE E^^^ du; ^ < . (75) 

It should be clear that even for At as small as l/m„ the integral is of order, , (we 
are extrapolating a non-relativistic formula to the relativistic region, which however 
should be allowed for an order of magnitude estimate), hence the quantity in the 
left hand side is at most of order, G\ ^ 10~^° , which means that the inequahty 
is violated only at At unreahstically small. It is thus difficult to observe deviation 
from the exponential decay in the short-time region of the neutron decay. 

The difficulty of observing deviation from the exponential decay in realistic par- 
ticle decay lies in a slow rise of the spectral function a{u;) near the threshold; 
a(a;) oc (a; — a;c)" with a > 2. 
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An implication of this consideration to the nucleon decay is that underground 
experiments performed so far can provide meaningful limits by using the limiting 
form 1 — Ft of the exponential law. This holds if the resolution time At > l/mjsf. 
A caveat for this conclusion is environment effect of hadronic matter, which will be 

discussed elsewhere. 

Finally, we mention that the infinite time resolution limit. At — 0, can be dealt 
with in a conventional approach of relativistic field theory assuming existence of 
the asymptotic field. The short-time limit Mt ^ of the formalism in Appendix 
however gives a peculiar behavior of non-decay amplitude if pcrturbative formula for 
the proper self-energy is used with ordinary rcnormalization. It is not clear how non- 
perturbative effect modifies this result. On the other hand, both the exponential and 
the late power behavior in the approach of Appendix is practically identical to the 
one in this section, if the theory is renormalizable. A great advantage of introducing a 
finite time resolution here is that one can discuss the non-renormalizable interaction 
on the same footing as the renormalizable one. 

A physical reason why the quadratic decrease in the short-time limit is obtained 
with a coarse time resolution is that a finite time resolution effectively cuts off the 
high frequency component, thus its outcome is similar to the behavior of a system 
of finite number of degrees of freedom. A field theory of infinitely many degrees of 
freedom may exhibit a behavior quite different from quantum mechanics of a finite 
system. Unfortunately, the true limit of t ^ e'^'^^/^^/M is difficult to reach in per- 
turbation theory. One may only probe the region oit ^ e^^'^^/^^/M in perturbation 
theory. 

6. S-wave decay 

We have discussed so far the case of Q not too small. The small Q limit may be of 
some interest in relation to the S-wave decay. The decay rate F goes with small Q as 
Q'+i/2 ^ where I is the angular momentum involved in the two particle decay. Thus 
the S-wave boson-pair decay (superrenormalizable model here) gives T/Q oc Q~^^'^ 
for a small Q and may become large. More precisely, the boson-pair decay in the 
Q — > limit gives 

r.^^0i^. (76) 

On the other hand, the P-wave fermion-pair decay (renormalizable model) gives 
r/Q oc Q^/^ . Thus, the pole locations move with the Q- value, quite differently 
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in the two cases of the S-wave and the P-wave decay, as shown in Fig. 6. For the 
exponential decay to exist, the imaginary part of the pole ~ V /2 <^ Q, and this 
condition may be violated for the S-wave decay of a very small Q <T value. 

In Fig. 7 we show an example of the behavior of time evolution for the S-wave 
boson-pair decay. The time evolution shown here corresponds to the parameter set 
marked in Fig. 6 by the crossed box. With a very small Q, the exponential form is not 
supported in any period. Indeed, the intermediate-time behavior in this case joins 
smoothly to a refined late-time formula of the power behavior. A refined calculation 
of the magnitude of the power law is necessary, because the Q limit gives a 
constant term in the formula of F{uj)~^ and the Q in eg . ([1^ is replaced by another 
linear term, although this constant is small with a small coupling. The late time 
formula is realized in this example at a time, a few times l/F, where F is given 
either by a naive perturbative formula or by a half of the imaginary part of the pole 
(both are roughly of the same order). It is however difficult to give a closed form 
formula of the transition time to the late power behavior, since the amplitude is not 
readily expressible in an analytic form. Nevertheless, the power of time in the late 
time formula is unambiguously determined by the threshold rise of the spectrum, 
(t{uj) oc (cu — cUc)" . It is an interesting question whether this type of behavior may 
be found in nature. 

In summary, we have formulated a non-perturbative definition of renormalization 
for the non-decay amplitude of unstable particle decay. We then discussed how the 
decay proceeds when viewed with finite time resolution. The exponential law usually 
holds as early as at ~ the resolution time. Exceptions may occur when the S-wave 
decay rate > Q, for which the power law may appear without the intermediate stage 
of the exponential decay. 
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Appendix. Non-decay amplitude and Green's function 

We first argue on general grounds that the non-decay amphtude should be finite in 
renormalizable theories. For this purpose let us recall the Lehmann type of spectral 
representation vahd in any models. Using a complete set of eigenstates of the total 

one may write for an arbitrary initial state |1), 

(l|e-*^*|l) = / duj p{uj)e-'^\ (77) 
p{u^)^Y.^i^-En)\{n\l)\\ (78) 

n 

Here |1) is not an eigenstate of the total H, while \n) is. Using the analytic function, 




that has the branch cut starting from z = 00^, the non-decay amplitude is given by 

(l|e-^^*|l) = — / du A(u; + iO+)e-'''K (80) 
2m J-00 

Thus, there exist two Fourier decompositions for the same object, this one and the 
one in the text. These two must be equivalent; F{u; + iO^) — A(a; -|- iO"*") . This 
further indicates the equality of the spectral functions; p{u;) — H{uj). We shall show 
how this spectral function is related to the usual Lehmann spectral function pl{u!^). 

If one takes the initial state to be an eigenstate of unperturbative Hamiltonian 
Hq (this automatically holds in the projector formalism), then one has a well-known 
relation to the time evolution operator; 

(l|e-'^*|l) = e-^^i*(l|C/(i,0)|l) = e-'^i*c(i), C/(i , 0) = e*^°* e"'^* . (81) 

As a physical requirement we impose that an unstable particle state |1) has no 
overlap with the true vacuum. This suggests to take as the initial state one particle 
state of the form, 

|1) = - i I d^xfo{x,ti) do ^{x,U) |0) , (82) 

since (0|</7(x)|0) = 0. Here foix ,ti) is the mode function appropriate for the initial 
state, (f{x ,ti) being the Heisenberg field operator, and |0) is the true vacuum of the 
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total Hamiltonian. The non-decay amplitude is thus written in terms of the Green's 
function, {0\(p{t) (p{ti)\0) ; 

c(() = e'"<"'> ^^^^tM J I rf3^. {0\TMxMx!m , (83) 

Ju>^ Z — u 

The Lehmann representation for the renormalized Green's function leads to 

^l\e-mt-U)\l) ^ 1 r rf^e--(*-*^)A(a;2 + iO+). (85) 

zni 2M J-oo 

If the spectral function is well convergent at a; — > oo, 

c{t) = e^^(*-*^) r duj ^^^^^^ PL (a;') e-^-(*-*>) . (86) 
The two spectral functions are then related by 

p(-) = ^^P.(c.^). (87) 

Presence of the factor {uj + M)^/2M in the last form of the non-decay amplitude 
casts a doubt on validity of exchanging the order of derivative and integration. There 
are however cases in which it is possible to introduce a regulator and deform the 
contour of uo integration such that the resultant renormalized formula is convergent 
even after taking the infinite regulator mass limit. 

The rest of arguments depends on existence of the asymptotic field corresponding 
to the unstable particle and perturbative expansion of the proper self-energy G(a;^), 

A(a;')-^ = a;' - - G{uj^) . (88) 

The proper self-energy G{uj'^) can be made finite by using the renormalized pertur- 
bation scheme (with necessary couter terms included as interaction terms). In effect, 
this renormalization is equivalent to changing the coefficient of uj"^ term in A(a;^)~^, 
to be later fixed by the normalization condition of the renormalized initial amplitude, 
Cii(O) = 1 , and at the same time using a twice-subtracted renormalized Gr{uj'^) 
instead of G{uj'^)] 

c«(t) = -i, r J r. .... . (89) 



2m 2M J-oo A^u^ -M^ -Gh(w2) ' 

G,(.') = 2.' ,3,f (90) 
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with M here the renormahzed mass. The spectral function a here coincides with 
(j{uj) of the projector formahsm to lowest order of perturbation. In lowest order of 
perturbation one may use cri{uj) already given for various models in the text. 

We note that eq. ( pUD is not well defined as an integral on the real axis, if the time 
derivative operation is performed within the u integral, due to a bad high energy 
behavior even in renormalizable models; (7(00) (x uj"^ as ^ 00. The perturbative 
infinity mentioned in is related to this ill defined formula (actually its expansion 
in terms of the proper self-energy). 

For the renormalized amplitude CR{t) it is more convenient to deform the integra- 
tion contour from the real axis. To any finite orders of perturbation there are only a 
finite number of branch cut singularities in the complex uo plane and deformation is 
possible passing through a finite number of Riemann sheets. Thus one may deform 
the contour into the lower half plane of S^tu > 0, where e"""^* gives a good convergence 
factor. In lowest order of perturbation one may use the complex integral formula 
encircling both the pole (Co) and Ci in Fig.l. Both contributions are well convergent 
even in the infinite regulator limit. The formula thus defined makes it posssible to 
estimate the short-time limit of the non-decay amplitude. 

We shall discuss the renormalizable model of fermion-pair decay. An explicit 
calculation in this model gives 

Stt UJ 

The Ci integral, in the t — > limit but for Mt ^ e"^'^^''^^, gives 



(l\e-'"'\l) ^ - iN e-'^^' r dye-y' = -iN — Ct:^ e"'^'* , (92) 
^ ' ' ^ IGtt^M Jo IGTT^Mt ' ^ ^ 

where is a constant to be determined by the wave function renormalization. The 
linear rise (oc 1/t) found here is surprising, but is inevitable if one can ignore the 
non-perturbative region of uj integral, uj > M e^'^'^l^'^ . The transition time from 
this very early to the exponential period may be estimated by equating the above 
amplitude to the exponential formula. The result is 
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This time scale is even much shorter than 1/M, and is beyond the reach of exper- 
imental search. A numerical evidence supporting the linear rise of the short-time 
limit is shown in Fig. 8. 

In the super-renormalizable model of boson-pair decay the short-time limit is 
different from the one in the renormalizable model above. It gives the linear decrease 
of 1 — 0\t] for the non-decay amplitude, however with a constant of the linear term 
different from F expected from the exponential law. 
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Figure Caption 



Fig.l 

Contours for integration of the non-decay amplitude. Dashed is the part of contour 
in the second Riemann sheet, with the dashed cross the location of the pole. 

Fig.2 

Comparison of large-time formulas. The solid curve is the real part of the non-decay 
amplitude obtained from an exact numerical integration along the contour Ci as in 
Fig.l, and the dashed one is the approximate formula in the text, the magnitude 
of amplitude computed in an arbitrary unit. A phase factor is taken out such that 
the real part of the exact amplitude becomes real and positive as t ^ oo. At times 
< 8.5 the exact real part becomes negative. 

Fig.3 

Three types of diagrams. 
Fig.4 

Non-decay probability for the fermion-pair model. In the inlet is shown an en- 
larged short-time behavior. Chosen parameters, both for this and figure 5, are 
g = I, Q 0.8977 Mo, M ~ 0.9977 Mq, with Mq the bare parent mass. A 
Lorentzian resolution function with Mq At = 20 is used. 

Fig.5 

Comparison of the short-time behavior with its limiting formula. The dashed curve 
is the pole approximation whose parameters are numerically obtained, while the 
dash-dotted one is the short-time limit in the text, eg (|6^) . 

Fig.6 

Pole location in the second Riemann sheet for the S-wave (crossed), and the P- 
wave(black circle) decays. Different Q values are computed for a fixed g = 1.0 and 
At = 20 in the unit of the bare mass Mq = 1. The dotted and the dashed curves 
are a fit of the T — Q relation in perturbation theory. The S-wave curve is better 
fitted with a formula, F = a^/Q~+~(], with constant a ,(3. The point marked by the 
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crossed box is the parameter set for Fig. 7. 



Fig.7 

S-wave non-decay probability. Chosen parameters for the boson-pair decay model 
are /i = 1.0 x Mq, M ~ 0.9988 Mq, Q ~ 1.07 x 10"^ Mq . The perturbative 
decay rate F ^ 1 x 10~^Mo. A Lorentzian resolution function with Mq At = 20 is 
used. 



Fig.8 

Short-time limit in the field theoretical approach outlined in Appendix. An exact 
numerical computation of Ci integral (the absolute value by the solid line and a 
dominant real amplitude by the dash-dotted fine) is compared to eq. ( P^ ) (the broken 
line) supporting 1/t rise in the fermion-pair decay model. Chosen parameters are 
the same as in Fig. 2. 
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